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A generalized continuous economic model is proposed for random markets. In this model, 
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that this type of markets reach their asymptotic equilibrium on the exponential wealth 
distribution. 
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1. Introduction 

In the last years, different techniques and models from statistical physics are being 
successfully applied to explain some real data observed in economy pQ. Concretely, 
it has been reported that western (capitalistic) economies can be divided into two 
distinct groups according to the incomes of people [5] . The 95% of the population, 
the middle and lower economic classes of society, arranges their incomes in an 
exponential wealth distribution. The incomes of the rest of the population, the 5% 
of individuals, fit a power law distribution. 
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Different models have been used to explain the origin of these wealth distribu- 
tions. From a macroscopic point of view, it is necessary to take into account that 
markets have an intrinsic random ingredient as a consequence of the interaction of 
an undetermined ensemble of agents which are performing an undetermined number 
of commercial transactions at each moment. A kind of models considering this un- 
knowledge associated to markets are the gas- like models [3]. These random models 
interpret economic exchanges of money between agents similarly to collisions in a 
gas where particles share their energy. In order to explain the two different types of 
statistical behavior before mentioned, different gas-like models have been proposed. 
On the one hand, the exponential distribution can be obtained by supposing a gas 
of agents that trade with money in binary collisions, or in first-neighbor interac- 
tions, and where the agents are selected in a random, deterministic or chaotic way 
|4|5|6j . On the other hand, the power law behavior can be simulated by introducing 
inhomogeneity in the gas by different methods, through the breaking of the pairing 
symmetry in the exchange rules [B] or through the saving propensity of the agents 

In this work, we consider a continuous version of one of such homogeneous gas- 
like models [S] which we generalize to a situation where the agents present a control 
parameter to decide the degree of interaction with the rest of economic agents. A 
zero value of the parameter will represent a null interaction of each agent with the 
rest of the ensemble and a value equal to unity will mean the total disponibility of 
each agent for the interaction with the whole gas. Under these circumstances, we 
are interested in explaining the appearance of the exponential (Gibbs) distribution 
for all the intermediary values, between zero and one, of the parameter value. This 
situation can represent the ubiquity of this distribution in many natural phenomena 
but in particular in the random markets. Thus, the rigorous analytical results here 
exposed will help us to enlighten the behavior of real data for income distributions 
but also to understand the computational statistical behavior found in such type 
of economic models which decay to the exponential wealth distribution in their 
asymptotic regime independently of the initial wealth distribution given to the 
system. 

2. The continuous gas-like model. Properties 

We consider an ensemble of economic agents trading with money in a random man- 
ner [4]. This is one of the simplest gas-like models, in which an initial amount of 
money is given to each agent, let us suppose the same to each one. Then, pairs of 
agents are randomly chosen and they exchange their money also in a random way. 
When the gas evolves under these conditions, the exponential distribution appears 
as the asymptotic wealth distribution. In this model, the microdynamics is conser- 
vative because the local interactions conserve the money. Hence, the macrodynamics 
is also conservative and the total amount of money is constant in time. 

The discrete version of this model is as follows (4] . The trading rules for each 
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interacting pair (rni,mj) of the ensemble of N economic agents can be written as 

m[ — € [nii + rrij), 

m'j = {I - e){mi + rrij), (2.1) 
i,j = l...N, 

where e is a random nmnber in the interval (0, 1). The agents are randomly 
chosen. Their initial money (mi, rrij), at time t, is transformed after the interaction 
in (toJ, m'j) at time t+1. The asymptotic distribution p/(m), obtained by numerical 
simulations, is the exponential (Boltzmann-Gibbs) distribution, 

Pf{m) ^ l3exp{-(3m), with f3 ^ 1/ < m >gas, (2.2) 

where pf{m)dm denotes the PDF (probability density function), i.e. the probability 
of finding an agent with money (or energy in a gas system) between m and m + dm. 
Evidently, this PDF is normalized, ||p/|| = pf{m)dm = 1. The mean value of 
the wealth, < m >gas, can be easily calculated directly from the gas by < to >gas = 

The continuous version of this model [8] considers the evolution of an initial 
wealth distribution pa{m) at each time step n under the action of an operator T. 
Thus, the system evolves from time n to time n + 1 to asymptotically reach the 
equilibrium distribution (to), i.e. 

lim T"(po(m)) -4p/(to). (2.3) 

n— >-oo 

In this particular case, pf{m) is the exponential distribution with the same average 
value < pf > than the initial one < pq >, due to the local and total richness 
conservation. 

The derivation of the operator T is as follows 8 . Suppose that p„ is the wealth 
distribution in the ensemble at time n. The probability to have a quantity of money 
X at time n + 1 will be the sum of the probabilities of all those pairs of agents 
(u,v) able to produce the quantity x after their interaction, that is, all the pairs 
verifying u + v > x. Thus, the probability that two of these agents with money (u, v) 
interact between them is Pn{u) *pn{v). Their exchange is totally random and then 
they can give rise with equal probability to any value x comprised in the interval 
(0, u + v). Therefore, the probability to obtain a particular x (with x < u + v) for 
the interacting pair (m, v) will be Pn(u) * pn{v)/(u + v). Then, T has the form of a 
nonlinear integral operator, 

p„+i(x) = Tp^ix) ^ [[ ^lil^fMiM audv . (2.4) 

J Ju+v>x u + v 

If we suppose T acting in the PDFs space, it has been proved [H] that T conserves 
the mean wealth of the system, < Tp >—< p >. li also conserves the norm (|| • ||), 
i.e. T maintains the total number of agents of the system, \ \Tp\ \ — \\p\\ — 1, that by 
extension implies the conservation of the total richness of the system. We have also 
shown that the exponential distribution Pf{x) with the right average value is the 



January 19, 2013 



4 Lopez-Ruiz, Shivanian, Abbasbandy, Lopez 

only steady state of T, i.e. Tpf = Pf. Computations also seem to suggest that other 
high period orbits do not exist. In consequence, it can be argued that the relation 
p.3p is true. We proceed to recall these properties and to establish some new ones. 

2.1. Properties of the operator T. Recall 

First, in order to set up the adequate mathematical framework, we provide the 
following definitions. 

Definition 2.1. We introduce the space of positive functions (wealth distribu- 
tions) in the interval [0,oo), 



Definition 2.2. We define the mean richness < x >y associated to a wealth dis- 
tribution y 6 Lf[0, oo) as the mean value of x for the distribution y. In the rest of 
the paper, we will represent it by < y >. Then, 



^0 

Definition 2.3. For a; > and y G L^[0, oo) the action of operator T on y is 
defined by 



with S{x) the region of the plane representing the pairs of agents (u, v) which can 
generate a richness x after their trading, i.e. 

S{x) = {{u, v), u,v > 0, u + V > x}. 

Now, we remind the following results recently presented in Ref. [3]. 

Theorem 2.1. For any y G L^[0,oo) we have that \\Ty\\ — (It means that 

the number of agents in the economic system is conserved in time). In particular, 
consider the subset of PDFs in _L^[0,oo), i.e. the unit sphere B — {y d L\^,oo), 
\\y\ \ = 1}. Observe that if y B then Ty G B. 

Theorem 2.2. The mean value < y > of a PDF y is conserved, that is < Ty >=< 
y > for any y G B. (It means that the mean wealth, and by extension the total 
richness, of the economic system are preserved in time). 

Theorem 2.3. Apart from y — 0, the one-parameter family of functions ya{x) = 
ae^"^, a > 0, are the unique fixed points of T in the space ^^^[0, oo). 



U{0}, \\y\\ < oo}. 



with norm 
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Proposition 2.1. The operator T is Lipschitz continuous in B with Lipschitz con- 
stant < 2. It means that if we take y,w £ B , then 

\\Ty^Tw\\<2\\y-w\\. 

Proposition 2.2. For any y e L^[0,oo) and m < n e N , it holds that 
(— l)'"(r"y)(™' G L'^[0,oo). It implies that if T^^y = y then y is a completely 
monotonia function. 



2.2. Other properties of T 

Evidently, our main interest resides in tlie action of the operator T in tlie subset of 
PDFs. We give here other properties for this restriction of T in B, although most 
of these properties are clearly valid for the whole space ^^"[0, oo). 

Proposition 2.3. Suppose that y & B then Ty is a decreasing function. 

Proof. Assume that 0:2 > and Xi < a:;2 then {(m, w), u,v > 0, u + v>X2}C 
{{u,v), u,v > 0, It + w > xi} and so 

^ [[ dudv > [[ dudv = Ty{.,) . 

Therefore Ty for every function y G -B is a decreasing function. This is a particular 
case of Proposition 12.21 for n — m — 1. □ 



Corollary 2.1. Suppose that y E B is a not decreasing function then Vn G N : 

T^'y + y. 



Proof. It follows directly from Proposition 12.31 □ 
Proposition 2.4. 

1. For every y,w E B we have \ \Ty — Tw\ \ < 2. 

2. For some members y,w G B, \\Ty —Tw\ \ — \ \y—w\\, hence T is not a contraction. 
Proof. 

1. It is clear that ||ry ~Tw|| < ||ry|| + \\Tw\ \ = 1 + 1 = 2. 

2. Take y = ae~"^ and w = fie^^^, with a, /3 > 0, which belong to B. In this case, 
we can easily see that \\Ty — Tw\ \ = \ \ae~"^ — I3e~'^^\\ = \ \y — w\\. □ 

Example 2.1. Take y{x) — (^i_^^yi and w{x) — which belong to B. By using 
Mathematica, it is seen that ||y-w|| = 0.407264 and \ \Ty-Tw\ \ = 0.505669. Then, 
\\Ty — Tw\ \ > \\y — w\\ which confirms again that T is not a contraction in B. 
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Conjecture: For any y G _B, we guess by simulation of many examples the 
truth of the following relation: 

{5e~''^ with 6 = 1/ <y>, 
or 
0^ when < y >= +oo . 

Let us observe that the above pointwise limit of T"y when n oo can be outside 
of B in the case that < y >= +oo. See the next example. 

Example 2.2. Take y{x) = (^i^r^^2 which belongs to B, with < y >= +oo. Evi- 
dently, T"y e B for all n. But it can be graphically guessed that lim„^oo T"-y{x) = 
0+ ^ B. See Figure [H 




Fi 



g. 1. (a) yix) = (b) Ty(x), (c) T^y{x). 



Theorem 2.4. Suppose that lim„_>oo WT'^uix) — fi{x)\\ = 0, and fi{x) is a continu- 
ous function, then fJ,{x) is the fixed point of the operator T for the initial condition 
y(x) G B. In other words, iJ,{x) = Se~^^ with 5 = 1/ < y > . 

Proof. First we prove that for each e > there exists M so that Vn > M: \\T"y — 
T/ij I < e. Since lim„^oo | \T"'y{x)—iJ.{x)\ \ = 0, we can say that 3N such that \/n > N: 
Ijjin-iy _^|| ^ g^2. Now, we choose M = N. Then, due to the Lipschitz continuity 
of T, we have 

Vn>M: ||T"y-r/i|| <2||T"-iy-/i|| <2.| -e. 

Now, by uniqueness of the limit, it implies that T/i = fj,. Therefore, for the initial 
condition y{x) G B, ^{x) is the fixed point of T, and by Theorem 12. 31 it means that 
^{x) = (5e~''^ with 5 = 1/ <y>. □ 

As particular examples, we present in Figs. [2] and [3] the graphical evidence of 
the convergence suggested in the former Theorem. Many other cases have been 
studied with a similar behavior. Remark that in Example 12.21 there is pointwise 
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convergence to the null function (fi — 0^) but there is no convergence in norm Li, 
in fact ||T"y - m|| = 1 for all n. 

Example 2.3. Take the Gamma distribution y{x) — xe~^, so that y £ B and 



^, then in this case ^{x) = 2^. We find numerically that 



W - Ml 



0.368226, \\Ty- ^JL\\ = 0.185608, \\T'^y - = 0.103225, \\T^y - = 0.061195, 
|r^y — ^11 = 0.037675, and so on. It is shown in Fig. [2] Then we can guess that 

lim„^oo \ \T'"^y - \ =0. 




(a) (b) 
Fig. 2. Plot of y{x) = xe~^ , T-iterates of y and fJ.{x) = i 

lir^y-Mll. 




e-2-. (a) II?; -Mil, (b) \\Ty - ^^\\, (c) 



Example 2.4. Assume now the rectangular distribution: y{x) = ^ if 2 < a; < 4, 
and y{x) = otherwise. So, y £ B and 5 — ^, then in this case fj,{x) — ie~3^. We 
find numerically that \\y - nW > \\Ty - > \\T^y- nW > \\T^y- and so on. It 
is shown in Fig. |3l Then we can also guess in this case that lim„_>oo | \T^y — m| | = 0. 




(a) (b) (c) 



Fig. 3. Plot of y{x) = i if 2 < x < 4, and y{x) = otherwise, T-iterates of y and tJ,{x) = ^e 3^. 
(a) II?/ -Mil, (b) ||T?;-m||, (c) ||T2?;-m||. 

Summarizing, the dynamical evolution of T in Lf[0, 00) is roughly given by the 
property: \\Ty\\ — , which divides the space L^[0,oo) in three regions: the 

interior of the unit ball B^ — {y £ L^[0,oo), \\y\\ < 1}, the unit sphere B, and the 
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exterior of the unit ball B^x = {y G L^[0,oo), \\y\\ > 1}. Accordingly, the recursion 

yn{x) ^T'^yoix), with yo{x) e L^[0,oo), 
presents four different behaviors depending on t/q: 

• If 2/0 e Bin, then lim„^oo ||2/ri(a;)|| = 0. 

• liyo & B and < yo > is finite, then lini„_i.oo \\yn{x) — Se^^^\\ — 0, with S = ^y^y ■ 
Observe that the succession ?/„ remains in B \fn. 

• li yo G B and < yQ > is infinite, then lini„_^oo ||2/n(2;) ~0+|| — 1. Observe that in 
this case the succession y„ remains in B \fn and the series yn does not converge 
in the norm || • || to its pointwise limit 0+ ^ B. 

• If 2/0 e Bex, then lim„^oo ||2/ri(a;)|| = oo. 

Hence, the only fixed points of the system are y = and Se~^^. They are asymp- 
totically reached depending on the initial average value < yo >, which determines 
the final equilibrium. We proceed now to show that this behavior is essentially 
maintained in the extension of this model for more general random markets. 

3. Generalized continuous model for random markets 

Let us observe that the action of operator T can receive a macroscopic interpretation 
respect to Eqs. (|2.1|) in the sense that, each iteration of the operator T means 
that many interactions, of order N/2, have taken place between different pairs 
of agents. As n indicates the time evolution of T, we can roughly assume that 
t « N *n/2, where t follows the microscopic evolution of the individual transactions 
(or collisions) between the agents (this alternative microscopic interpretation can 
be seen in [10]). 

Let us think now that many of the economical transactions planned in markets 
are not successful and they are finally frustrated. It means that markets are not 
totally effective. We can refiect this fact in our models in a qualitative way by defin- 
ing a parameter A G [0, 1] which indicates the degree of effectiveness of the random 
market. When A = 1 the market will have total effectiveness and all the operations 
will be performed under the action of the random rules (j2.ip . The evolution of the 
system in this case is given by the operator T. When A — 0, all the operations 
become frustrated, there is no exchange of money between the agents and then 
the market stays frozen in its original state. The operator representing this type 
of dynamics is just the identity operator. Therefore, we can establish a generalized 
continuous economic model whose evolution in the PDFs space is determined by 
the operator T\, which depends on the parameter A as follows: 

Definition 3.1. 

Txy{x)^{l~X)y{x) + \Ty{x), with Ae[0,l]. 

Observe that the parameter (1 — A) can also be interpreted as a kind of saving 
propensity of the agents, in such a way that for A = 1 they do not save anything and 
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they game all their resources, and for A = they save the totality of their money 
and then all the transactions are frustrated and the market stays in a frozen state. 

3.1. Properties of the operator Tx 

Now, we proceed to present the properties of the operator T\, which shows a dy- 
namical behavior essentially similar to the behavior of T. Concretely, one of our 
main results is that the exponential distribution is also the asymptotic wealth dis- 
tribution reached by the system governed by T\, independently of the effectiveness 
A of the random market. 

Let us observe that Tx = ^ for A = and T\ — T ioi X — 1, where / is the 
identity operator. 

Proposition 3.1. T\ conserves the norm, i.e., for each y ^ B, we have T\y G B. 

Proof. Consider y{x) > 0, then T\y{x) > 0. The norm is: 

\\nyix)\\ = (1 - \)Mx)\\ + X\\Tyix)\\ = (1 - A) + A = 1. □ 

Proposition 3.2. T\ conserves the average value ofyGB, i.e. < y >=< T\y >, 

where < y > represents the mean value < x >y as expressed in Definition \2.2\ 

Proof. Take y £ B, then 

< Txy >= {I - X) < y > +X < Ty >= {1 - X) < y > +X < y >=< y> . □ 

Proposition 3.3. The operator Tx is Lipschitz continuous in B with a Lipschitz 
factor a such that a < 1 + A. 

Proof. Suppose that y,w E B, then 

\\Txy-Txw\\ = ~ X)y + XTy ^ {1 - X)w - XTw\\ < 

< (1 - X)\\y -w\\+ X\\Ty - Tw\\ < 

< (1 - A)||y - w\\ + 2X\\y - w\ \ ^ 

< (1 + A)||y-^«||. □ 

Theorem 3.1. For any X G (0,1), the operators T and Tx have the same fixed 
points. 

Proof, (i) Suppose y is a fixed point of T, i.e. Ty = y. Then we have Txy = y, 
because Txy = (1 — X)y + XTy = y — Xy + Xy = y. 

(ii) Now suppose y is a fixed point of Tx, i.e. Txy = y. Then we have Ty = y, 
because XTy = Txy - {I - X)y ^ y - y + Xy ^ Xy. □ 
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Corollary 3.1. The function y{x) = and the family of exponential distributions 
ys{x) — Se^^^, 6 > 0, are the only fixed points ofT\ in ij'"[0, oo), with A G (0, 1]. 

Proof. It follows from Theorems 13.11 and 12.31 □ 

Theorem 3.2. Suppose that for a given A G (0,1) we have lim„_j.oo \ \T^\y{x) — 
^(a;)|| = 0, with iJi{x) a continuous function, then fJ.{x) should be the fixed point of 
the operator T\ for the initial condition y{x) E B. In other words, — Se~^^ 

with S = 

<v> 

Proof. Identical to the proof of Theorem 12.41 by changing T by T\ and taking into 
account Theorem 13. II □ 

Example 3.1. Take again the Gamma distribution y{x) = xe~^ , so that y € B and 
S — ^, then in this case iJ,{x) = ^e^2^. For A = 0.5, we find numerically that \\y — 
H\\ = 0.368226, \\Txy-ii\\ = 0.273011, \\T^y~ti\\ = 0.206554, \\T^y-fi\\ = 0.158701, 
and so on. It is shown in Fig. |4l Then we can guess that lim„^oc WT^y ^ t^W =0. 
Observe that in this case Txy can be a non-decreasing function. 




4. Conclusions 

In this work, a continuous economic model recently introduced [8] has been general- 
ized. This model takes into account idealistic characteristics of the markets, where 
agents interact by pairs and exchange their money in a random way. Also, the model 
implements a parameter that gives an idea of the effectiveness of the agents when 
trading between them. A perfect effectiveness means the total availability of the 
agents wealth to be gamed in their transactions and a null effectiveness gives rise 
to a frozen market where each agent keeps intact his money. We have shown in a 
rigorous manner that it does not matter the degree of effectiveness of the market 
for the final statistical result. In all the cases, these random markets evolve toward 
their asymptotic equilibrium, that is, the exponential wealth distribution. 
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